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Spectral weight functions are easily obtained from two-point correlation functions and they might 
be used to distinguish single-particle from multi-particle states in a finite- volume lattice calculation, 
a problem crucial for many lattice QCD simulations. In previous studies, it is shown that the 
spectral weight function for a broad resonance shares the typical volume dependence of a two- 
particle scattering state i.e. proportional to 1/L in a large cubic box of size L while the narrow 
resonance case requires further investigation. In this paper, a generalized formula is found for 
the spectral weight function which incorporates both narrow and broad resonance cases. Within 
Liischer's formalism, it is shown that the volume dependence of the spectral weight function exhibits 
a single-particle behavior for a extremely narrow resonance and a two-particle behavior for a broad 
resonance. The corresponding formulas for both Af and T[' channels are derived. The potential 
application of these formulas in the extraction of resonance parameters are also discussed. 
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I. INTRODUCTION 



Lattice Chromodynamics (lattice QCD) has provided 
us with a genuine non-perturbative theoretical frame- 
work to study low-energy physics in strong interactions 
from first principles. Combined with numerical Monte 
Carlo simulations, lattice QCD has become a crucial 
player in hadron spectroscopy and other studies in the 
field. 

In a typical lattice simulation, one computes a two- 
point function of an appropriate operator with a defi- 
nite quantum number via Monte Carlo simulation and 
the energy eigenvalues in that sector are then extracted 
from the time-dependence of the correlation function. 
In fact, the correlation function of an operator O may 
be written schematically as: C{t) = (C(t)C(O)) = 
^„ VF„(i?„)e~-^"* where E^s are the energy eigenval- 
ues of the QCD Hamiltonian for the corresponding state 
n with the quantum numbers associated with the opera- 
tor O. The functions W„(£'„) are the so-called spectral 
weight functions for the state n. Since lattice simulations 
are all performed in a finite volume, the eigenvalues En 
are all discrete in nature. This poses a question on the 
particle nature of the corresponding state. If in the infi- 
nite volume, single-particle and multi-particle states arc 
different in the sense that the spectra for single parti- 
cle states are discrete while that for multi-particle states 
are continuous starting from the multi-particle thresh- 
old. In a finite volume, however, it is not an easy task 
to tell whether a state with a particular energy En is a 
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single- or a multi-particle state since all states have be- 
come discrete. This distinction is not only for academic 
purposes but also of phcnomenological importance as well 
if one recalls that quite a number of new resonance-like 
structures have been identified rather close to the multi- 
particle threshold of known hadrons. To cope with this 
problem, it has been suggested that the volume depen- 
dence of the spectral weight functions Wn{En, L), which 
is also a function of the box size L, may be used as a 
probe to distinguish single- and multi-particle states [l[ . 
With an appropriate normalization, it is argued that the 
spectral weight function for a single-particle state has lit- 
tle volume dependence while that for a two-particle state 
will exhibit a 1/L^ behavior. 



In our previous study [J , it was shown that the above 
mentioned criteria becomes invalid for a broad resonance. 
Instead of showing little volume dependence, the spectral 
weight function for a broad resonance will exhibit a 1/L^ 
behavior, typical for two-particle states. However, our 
previous study does not apply for extremely narrow res- 
onances. In this paper, we re-examine the problem within 
Liischer's formalism (see Ref. Q), with an emphasis on 
the narrow resonance limit case. A generalized formula 
is found which incorporates both narrow and broad reso- 
nance cases. The conclusion we reached is exactly analo- 
gous to what we discovered in a model study Q , namely 
that for a broad resonance, the spectral weight behaves 
like a multi-particle state while for an extremely narrow 
resonance, a single-particle behavior realizes. We also 
point out the possibility to extract the resonance param- 
eters from the spectral weight function. We closely follow 
the discussion of our previous work [2] . The reader is re- 
ferred to that reference for notations and further details. 



II. THE MODEL IN THE INFINITE VOLUME 
AND THE SIGNATURE OF A RESONANCE 

In the infinite volume, consider a non-relativistic quan- 
tum mechanical model with Hamiltonian given by: 



H = — ^—V^ + V{r 
2m 



(1) 



where the potential ^(r) is zero for r > a with some 
a > 0. Wc now discuss the energy eigenstates satisfy- 
ing: H'^{r) = E'i^{r). One can expand the eigenfunction 
as: *(r) = tpim{r)Yim{n) with: r = rn and ■0/„i(r) is 
the radial wave- function satisfying the radial Schrodinger 
equation: 



dr'^ 



2 d l{l + l) 
r dr r^ 



+ k' - 2mVir) ] ^Pi„,{r) = 



(2) 

where E = k'^/{2m) being the energy eigenvalue of 
the state. There exists a unique solution to the ra- 
dial Schrodinger equation that is bounded near the ori- 
gin which is denoted as: ui(r]k) and the general solu- 
tion to the radial Schrodinger equation has the form: 
'4'imif) = bimUi{r; k) with some constant 6;,„ to be fixed 
by other conditions. 

In the region r > a where the interaction vanishes, the 
solution ui (r: k) are expanded in terms of spherical Bessel 
functions: 



ui(r:i 
^:0 



ui{r; k) = ai{k)ji{kr) + jii{k)ni{kr) 



(3) 



The coefficients ai{k) and I3i{k) have simple relation with 
the scattering phase shift: 



^2.M^) = ^'(f) + ^/3'(fc) ,ot5,{k) = 



ai{k) 



(4) 



ai{k)^ipi{k) ' '' ' pi{k) 

In the low-energy limit /c — > 0, one normally defines: 

a? = lim k^ai{k) , /S^* = lim k-^-^f3i{k) , (5) 



fc^O 



fc->0 



and the threshold parameters ai = (3f /a'}. In particular, 
oo for Z = is referred to as the s-wave scattering length. 
The threshold parameters ai are important because they 
characterize most behaviors in low-energy scattering pro- 
cesses. For example, we have: 



Si{k) ~ a,fc2'+i + 0(fc2'+3) , (mod tt) 



(6) 



To fix the normalization of ui{r; fc), we impose the con- 
dition: (si 



ui(r; fc) ~ — sin(fcr — lTr/2 + SA 
kr 



(7) 



Another way of writing the asymptotic behavior of the 
radial wave function in Eq. ([3]) is as follows: 



1 



uiir-k) c^ -[Ai[k)e'''''- + Bi[k)e 



ikr~\ 



(8) 



where the functions Ai(k) and Bi(k) are related to the 
coefficients ai{k) and Pi{k) via 

Mk) = {-i)\lii-iai)/{2k),Bi{k) = {+i)\Pi+iai)/{2k). 

(9) 
When regarded as functions of the energy E = fc^/(2m), 
we will write them as Ai{E) and Bi{E), respectively. 
Note that for real positive energy, we have: 



Ai{E)^Bi{Ey , E>0,Ee 



(10) 



The condition for a narrow resonance is that, the coef- 
ficient Bi (E) has a complex zero on the second Riemann 
sheet ai E = Ei, ~ Jr/2, close to the real axis: 



(11) 



Bi[E,-i-]=Q, 



with r being the physical width of the resonance. In this 
case, the scattering properties of the system is mainly 
governed by the position of this resonance pole, if the 
scattering energy E is very close to the resonance energy 
E^,. In fact, since the function Bi(E) vanishes at _E* — 
ir/2, for real energy E close to i?*, we may approximate 
the function by: 



Bi{E)c^hi{E-E^+iV/2), 



(12) 



for some complex constant hi. Therefore, for positive real 
energy E close to i?*, the radial wavefunction looks like 

ui{r- k)~- [b*i{E -E^- iVl2Y^''' + h.c\ . (13) 

This then shows that the phase-shift near the resonance 
point is given by 



.2..,(i.)„,2..;-7 g-^^-»r/2 

E-E^+ ir/2 



{*) 



where the quantity 6, is given by 



e'^'' =i-iy+'{b:/bi) 



l+ll 



(14) 



(15) 



This is a constant phase (meaning that it does not de- 
pend on the energy substantially in the resonance region) 
which depends only on the position of the resonance. The 
energy dependent part of the phase shift is given solely 
by the factor in the big parenthesis in Eq. (|14p . It is 
seen that, for extremely narrow resonance, the phase shift 
jumps by tt when the energy passes through E^, within 
an energy range of order of several F. 

Another important point has become clear from the ex- 
pansion in Eq. P^ . If we set the energy E to be exactly 
the resonance energy _E^, the expression ui is vanishing 
if F becomes infinitesimally small. This means that, to 
properly normalize the wave-function, therefore, the co- 
efficients bi has to be of order 1/F. 



III. THE MODEL ON A TORUS 

Now we enclose the system we discussed in the pre- 
vious section in a large cubic box of size L and impose 
the periodic boundary condition in all three spatial direc- 
tions. The potential itself is also periodically continued 
to Vl{y) = Enez3 ^(|r + nL|). We define the so-called 
"outer region" as 



17 = {r : |r + nL| > a, for aU n e Z^} . 



(16) 



This is the region where the potential vanishes identically. 
We assume L » 2a so that the outer region admits free 
spherical wave solutions (asymptotic states). Following 
the discussion in our previous paper, we may write the 
true eigenfunction as: 



*(r; fc) = ^ bunui{r- k)Yi^{n) . 



(17) 



where the coefficients bim are to be determined by bound- 
ary conditions and normalization. 

In the outer region, the eigenfunctions can also be ex- 
panded in terms of singular periodic solutions (SPS) of 
Helmholtz equation. These are defined as 



Gi™(r; k) = yi^{V)G{v- k) , G(r; k) ^ ^Y. 



7,3 A^ p2 „ fc2 ■ 

p 

(18) 
with 3^;m(r) = r^yimi^r) being the harmonic polynomi- 
als. We may write 



*(r;A:)|r6n = ^ Wi™Gim(r; fc^) 



In the meantime, the outer solution can also be ex- 
panded in terms of spherical harmonics and the spherical 
Bessel functions ji{kr), ni{kr) and the matrix elements 
■Mi'm'-.imikf). The explicit expressions can be found in 
Ref. 13. 

It is also convenient to introduce the cubic version for 
the functions Gim defined by Eq. ([T8|) : 



C(r;fc)=3^^?(V)G(r;fc) 



(20) 



;(n, 



(19) 



where yi\ (r) is the cubic harmonics for a particular rep- 
resention F of the cubic group. The explicit expressions 
are listed in Appendix \X\ at the end of this paper. 



A. Liischer's formula in Af and Tj channels 

In this subsection, we will derive Liischer's formula in 
the A^ and T^ channel (corresponding to s-wave and 
p-wave in the continuum, respectively). Using similar ar- 
guments and the cubic harmonics given in the appendix, 
it is easy to generalize the results to other channels. 

We first examine the A^ sector which is the analogue of 
s-wave scattering. A good approximation for the s-wave 
dominated eigenfunction can be written as a superposi- 
tion of / = and I =^ 4 cubic harmonics with the s-wave 
component much larger than that of g-wavc. So, we may 
write the eigen-function in Af sector as: 



*(^^)(r;fc) =6o«o(r;fc)yJ^^'(r!r) + M4(r;fc)ri''^^(f7r) + --- , 

I 



(21) 



with 1 64 1 <C \bo\ in the large volume limit. In the outer we have 
region where the interaction vanishes, using relation ^ 



¥^i'>{r;k)\ren = 6o[«oJo(fcO + /?ono(fcO]>^o^^^(^r) + &4[a4j4(fcr) + ^4n4(fcr)]yi^^^(fir) 



(22) 



r 



On the other hand, in the outer region f2, we may Helmholtz equation defined in Eq. (|20l) : 
also expand the wavefunction in terms of cubic SPS of 



J 



^i^r){r;k)Un={^]v„ 



G^.'^\r;k)+'^A'^\r;k) 



(23) 



Matching the above two expansions in the outer region 
and using the exphcit expressions for the cubic harmon- 
ics, one arrives at the following result (i.e. Liischer's 
formula) for the Af sector: 



(cot(5^°'' — TOoo)(cot(5^ ' — 17144) = ?Ti4omo4 



(24) 



For the sector T^ , which is the cubic analogue of the 
vector channel, we have a similar expression: 



vff i"^(r; k) = &iui(r; k)Yl^^^\n,) + b^usir; k)Y^^/\nr) + ■■■ , 



(25) 



r 



where we have only written out the third component of of cubic SPS of Helniholtz equation: 
the vector wavefunction. We may also expand in terms 



J 



* 



(Tf) 



(r;fc)|,. 



en 






^i;3 (i";fc) + ^S3;3 (r;^) 



(26) 



Using the explicit expressions for the cubic harmonics 
and cubic SPS to match the above two equations we ob- 
tain Liischer's formula in the T^ sector: 



{cot 5^^' — mii)(cot j'^^ — m-i^) = 7711377131 



(27) 



Since numerically we have: 771 n = 77100, therefore if 
higher momentum contaminations (i.e. I > 4 for A~[ 
and I > 3 for T^) are neglected, the formula in the T^ 
channel has the same form as the one for Ai channel. 



B. Spectral weight function in A'{ sector 

As we pointed out in our previous study Q, the spec- 
tral weight function has to be normalized in a proper 
manner so that the single-particle states has a conven- 
tional behavior. For an interpolating operator 0{t), the 
state O^(0)|0) will generally be a superposition of a tower 
of energy eigenstates: {|i?)} with energy eigenvalue E. 
The spectral weight for a particular energy state \E) ob- 
tained from the two-point correlation function of 0{t) is 
given by: 



W{E) = \OiE)\'^\{E\O\0)m' 



(28) 



where 0{E) = {E\0^ {0)\0) is the overlap of the exact 
energy eigenstate \E) with the state |C''(0)|0). 

The exact energy cigenfunction (ri, r2\E) only depends 
on the relative coordinate r = r2— ri and does not depend 
on the center-of-mass of the two particles. Therefore, in 
terms of the wavcfunctions associated with the relative 



motion, the normalization conditions for the A^ and T^ 
channels become: 



d3r|*('4+)(j..^.)|5 



d^r 



Ts 



T3 
'^f''\r;k) 



1 



(29) 



*r"\r;fc) = -1^,:,, (30) 



where vl/('^i ' (r; k) and ^^ (r; k) is the wavefunction for 

the A^ and Tj" channel respectively. 

Let us first investigate the normalization condition in 
the A'^ channel. The integral over the torus in Eq. ((29)) 
runs over two regions: the inner region B = {r : r < 
a, mod L} where the potential is present and the outer 
region J7 where it vanishes: 



rf3j.|^(^f)(r;fc)|2_j_ / d3r|*('^^)(r;fc)|2 



1 

1? 



(31) 



If we assume that the wavefunction in this sector is dom- 
inated by the s-wave component such that 



vI/(^f)(r;fc)|r6B 
¥^^\r;k)\ren 



bouoir] k)Yoo{n) , 

47r\ 

-77 I voGoo{r;k) , 



(32) 



Matching the above two wavcfunctions at r = a will fix 
the ratio bo/vo- Therefore, wc may write: |6op = CI'^^oP 
with C being some positive constant with the dimension 
of k^^. Therefore, we may express the normalization 
condition as 



1^0 I 



C / r'\uoir;k)\'dr^ 





47r 
k 



r 



— / d'r\Goo{r;k)\' 



1 

l3 



(33) 



If wc drop the first term in the bracket in the above equa- 
tion and approximate the second integral by an integral 
over the whole torus, we obtain exactly the same for- 
mula as we derived in our previous study (see Eq. (38) in 
Ref. [2| and the discussions about it). Normally, if there 
were no sharp resonance in the scattering, neglecting of 
the first term is justified since the second term is pro- 
portional L^, much larger than any normal term in the 
large volume limit. However, when an extremely narrow 
resonance is present, this is not the case, as we remarked 
at the end of Sec. HT] 

The first term in the normalization condition involves 
the exact radial wavefunction in the interaction region 
which seems to depend on the details of the potential. 
However, although generally speaking we do not know 
the exact form of the inner wavefunction uo(f;fc), the 
integral /„ r'^\uQ{r;k)\'^dr can be related to the corre- 
sponding phase shift via an exact formula: [9| 



r |uo(r; fc)|' dr 



dSojk) 
dk 



k 



sin2(A:a + ^o) ■ 
(34) 



Now, if there exists an extremely narrow resonance at a 
particular energy Ei, = k1/{2m), then close to the re- 
soance we have cotJo — —{E — i?*)/(r/2). This implies 
that: 



d5o{k) 
dk 



k dSo 
m dE 



r/2 



{E^ - Ef + r2/4 



(35) 



It is seen that, close to the resonance energy, the deriva- 
tive d6a/dk is proportional to l/F in the narrow reso- 
nance limit (i.e. when E c:i Ei, and F ^- 0). Therefore, 
in the narrow resonance limit, such a contribution can 
become comparable or even larger than the second term 
that is proportional to L^ in a finite volume. All the 
other terms in the integral /„ r^ 1110(7"; k)\'^dr are bounded 
in both the narrow resonance and the large volume limit. 
Therefore, taking into account of this possibility, the nor- 
malization condition for the wavefunction is written as 



l^'ol 



cr 



{E,~Ef+T^/A 



An 



r 



d\\Goo{r;k)\' 



n 



1 



(36) 



where ^ and C > are two constants that are regular in 
both the large volume and the narrow resonance limit. 

The integral in Eq. (j36|) has been computed in our 
previous paper with the result [2|: 



A-TT 

d'r\Goo{r;k)\' = —F\k'), 



T3 



fc2 



with the function F{k'^) defined by 



F(fc2 



1 V fip') 

P 



l^'ol 



cr 



(37) 



(38) 



In the large volume limit, the derivative of this function 
has the approximate form 



F'{k^) ~ ^ cot 5o{k) + ^^ csc^ do{k) . (39) 



Collecting everything together, we obtain the following 
normalization condition in the A^ sector: 



iE^-E)^+V^/A 2A;3 



—-T cot Sa{k) + csc^ 6o{k) 



/cAp2 



1 

Z3 



(40) 



As discussed above, the second term in the bracket in the above formula is singular in the narrow resonance limit 



(0(l/r)) while the fourth term is singular in the large 
volume limit {0{L^)). All other terms are regular in both 
limits. 



Following similar steps as in our previous paper [31 , the 
speetral weight funetion in the A^ ehannel now takes the 
following form: 



J 



w^-^^Xe) 



47:\ipL{k')\yk 






1 

2P 



cot(5o(fc) 



fcAp2 



■Soik) 



(41) 



r 



where the function ipL^k'^) has been introduced in Ref. 0| 
which has little volume dependence. The main difference 
between this formula and the one obtained in our pre- 
vious paper is the appearance of the term (Cr)/[(£'* — 
E)'^ +r^/4] in the denominator. This term becomes sin- 
gular in the narrow resonance limit (proportional to l/F, 
as E = Ei,). If the resonance is not extremely narrow, 
this term is non-singular and can be absorbed into the 
constant term ^ and we recover our old result: the res- 
onance acquires a spectral weight function proportional 
to l/L^, typical for two-particle states. However, if the 
resonance is extremely narrow such that: 



C/r» 



1 



A:Ap2 



(42) 



the denominator is dominated by the term containing 
r and the whole spectral weight function has little vol- 
ume dependence, typical for a single-particle state. It is 



also noted that, in this limit, the spectral weight func- 
tion is proportional to F, the physical width of the res- 
onance. This is expected since when the resonance be- 
comes extremely narrow, the coupling of the resonance 
with the two-particle final states also becomes infinites- 
imally small (and is proportional to the width as it 
should). To summarize, Eq. (PT|) is a generalization of 
our previous formula in the sense that it incorporates 
both the large volume and the narrow resonance limit. 



C. Spectral weight function in T-^ sector 



Following a similar treatment as in the A[ channel, 
we can obtain the relevant formulas in Tj^ channel. We 
construct the state using an interpolating operator of the 
form: 



J 



$) =Ot(o)|o) = ^|.(P)P, [^t(p^o)^|(-P,0)-*t(p^o)^I(-P,0)l |0) 



(43) 



r 



where we have only taken one component (the third com- 
ponent) of the wavefunctions for a vector state. After 
factorizing out the center of mass motion, the normaliza- 
tion condition expressed in terms of the wavefunction for 
the relative motion becomes 



d^rl^fi 



(r;fc)|^ 



d3r|vl/fn(r;fc)p = -l 



(44) 



where the first integral is over the interaction region 
while the second is over the outer region where interac- 
tion vanishes identically. For the p-wave dominant wave- 



function, we have ^ 



(Tr 



(r;fc)|. 



5iWi(r;fc)Yio(rir) 



for the wavefunction in the interaction region and 



* 



{TD 



(r;fc)|reo — (47r/A:^) i;iGio(r; fc) in the outer re- 



gion. [l3| Matching the two wavefunction at the bound- 
ary yields the relation |6ip = C'l'^iP foi' some positive 
real constant (' > 0. The first integral in Eq. (|44|) can 
again be estimated in the same manner as that for Af 
sector: 



r |ui(r; fc)| dr 



dSijk) 
dk 



— sin2(A:a + (5i) . 
k 



Assuming that there exists a resonance aX E — E^, in this 
channel and for energies that are close to this resonance 
we have: cot S{E) ~ -{E - ^*)/(F/2). Following simi- 
lar derivations, we arrive at the following normalization 
condition for an energy eigenstate whose energy is in the 
resonance region: 



I^ll 



e' 



C'T 3 

{E^ - EY + r2/4 ^ 2fc3 



cot(5i(fc) 



fcAp- 



■csc^(5i(fc) 



1 



r 



(46) 



where ^' is some real parameter which is regular in both normalization condition in the T^ sector, the spectral 
the narrow resonance and large volume limit. With this weight function takes the following form: 



w^^^\e) 



lQn\a + k^^L{k^)Y/{ik^) 



^' + {E,-Ey^+ry4 + T^ ^°^ '^i (^^ + m^ 



esc 



Siik) 



(47) 



where the constant a is defined by 



(27r)3 



(fP^(P) , (48) 



Inspecting Eq. (jTZ]) we find that, the behavior in the 
T^ sector is the same as in the A^ sector qualitatively. 
The second term in the denominator is singular in the 
narrow resonance limit while the fourth term is singu- 
lar in the large volume limit. Which of these two terms 
dominates the spectral weight function depends on the 
size of the following two quantities in the denomina- 
tor: C'/r and l/{kAp'^). For an extremely narrow reso- 
nance: C'/r ^ l/(fc^P^) and the spectral weight func- 
tion is almost volume independent, a typical behavior 
for a single particle state. On the other hand, if we have 
C'/r < l/(fcAp2) the spectral weight exhibits a 1/L^ 
behavior, typical for a two-particle state. 



IV. CONCLUSIONS 

In this paper, the volume dependence of the spectral 
weight function is analyzed with the emphasis on the 
possible narrow resonances. Generalizing our previous 
results on this subject, we obtain a formula which in- 
corporates both the large volume limit and the narrow 
resonance limit. This is summarized in Eq. (j4ip for the 
A^ channel and in Eq. (|T7)) for the T^ channel. If so 
desired, formulas for other channels can also be derived 
similarly using the cubic functions listed in the Appendix. 
These formulae for the spectral weight functions exhibit 
the following general feature: near a resonance the spec- 
tral weight shows single-particle or mult i- particle behav- 
ior depending on whether the resonance is narrow or 
broad in the corresponding box. This is controlled by 
the size of the two quantities (/T and l/(fcAp2) that ap- 
pear in the spectral weight function. If the former one 
dominates, a single particle feature (an almost volume 
independent spectral weight function) realizes; if the lat- 
ter one dominates, a two-particle feature emerges with a 



l/i"^ dependence on the volume for the spectral weight 
function. Therefore, using the general Liischer's formal- 
ism, we have shown this general feature for the spectral 
weight function, just as we discovered it in a solvable 
model in Ref. [J]. 

Furthermore, the spectral weight function has a gen- 
eral dependence on the volume of the box and the phys- 
ical width of the resonance in the resonance region. As 
we pointed out in our previous study Q , this opens up a 
possibility to extract the resonance parameters from the 
spectral weight functions obtained in Monte Carlo simu- 
lations. It is readily verified from Eq. (|4T|) and Eq. (|47)) 
that, close to the energy of a resonance, both spectral 
weight functions satisfy the following functional form: 



1 



W{E,L) 



A{E) 



C{E) 

{E^-Ef+T'^/4. AE 



B{E)T 



(49) 



where A{E), B{E) and C{E) are smooth functions of the 
energy in the resonance region and they do not contain 
explicit volume dependence. The explicit volume depen- 
dence only comes from the term containing AE and the 
information about the resonance (i.e. Ei, and F) is en- 
coded in the second term of the above equation. If one 
performs simulations on a series of volumes with the ener- 
gies close to an underlying resonance and fits the spectral 
weight functions, it is in principle possible to extract the 
resonance parameters. Of course, the feasibility of this 
remains to be checked in realistic simulations. 



Appendix A: Cubic harmonics for various 
irreducible representations 

In this appendix, we list the cubic harmonics used in 
our calculation for various irreducible representations of 
the cubic group. We denote these cubic harmonics as 
3^; j (r) where F is the label for the irrep, I is the cor- 
responding angular momentum label (for the rotational 
group) and i — 1, ■ ■ • ,dim{T) denotes different basis 
functions of the irrep F. If F were one-dimensional, the 



label i will be dropped. The functions 3^;.j (r) defined 
here are homogeneous functions of order I . We will also 
use the notation: 



>l-'(^r)-r-'.4y(r) 



(Al) 



nary spherical harmonics. Note that these functions are 
orthogonal to each other but some of them are not prop- 
erly normalized. The table for these harmonics can be 
found in Ref. Q- 



to designates the corresponding angular functions. It is 
easy to see that they are linear combinations of ordi- 

I 

^ The A+ channel: 



yi''^\r) = y^o + y ^(3^4,4 + 3^4,-4) « L^ + y^ + z^- \A , 



(A2) 



«®° The A:; channel: 



,iA; 



y^^''{r)=xyz, 

y!,^^\r) = xyz (x'' + y^ + z^- ^A 



(A3) 



^ The E+ channel: 



h'^:\r 



)-{z'-\{^' + v')) ,yiPir)-{x'-y') 



y^^(r) 

y'ePir) 

yS\^) 



,2.^;(S+)/ 



■^' -\{^' + y') - \r^ ■ 3^^r^(r)) , 3^if ^(r) = [x^ - y' - \r^ ■ y^ ' ^(r) 



' ^' %^r--y'S\-)-Y--y'£\-) 



X -y 



(A4) 



The T{ channel: 



3^S'^(r)=3^io(r)cxz, 
yg^^iv)^yMr^)o,z(z'^\r' 



y^;^p{v) = z[z^ 



10/2 3 2\ 2 34 
z r \ r — —r 



7(^1 )' i^\ ^ J ^i , „.4 3 ^^2 , .,2^2 



yi;r(^)-A^' + v--A^' + ry 



(A5) 



^s- The T+ channel: 






yg^y(r) = xy(x' + y'-^-{x'+yr 



(A6) 



Note that for the irrep T^ and T^, the basis functions 
is formed by 3 functions. Here we only hst the ones that 
correspond to the z components. Other two sets can be 
obtained by cychcly swapping the coordinates. It is also 
noted that at Z = 5, we have two sets of basis (one of 
the basis function is denoted by a prime) since the irrep 



T-^ appears twice in the decomposition of / = 5. Similar 
situation occurs for irrep Tj at I = 6. Since the function 
G{v]k) is cubic invariant, it is easy to verify that Ql.^ 
defined in this Appendix indeed forms a basis for the 
corresponding representation T of the cubic group. 
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